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ABSTRACT 

A set of co-aligned high-resolution images from the Atmospheric Imaging Assembly (AIA) on board the Solar 
Dynamics Observatory is used to investigate propagating disturbances (PDs) in warm fan loops at the periphery of 
a non-flaring active region NOAA AR 11082. To measure PD speeds at multiple coronal temperatures, a new data 
analysis methodology is proposed enabling a quantitative description of subvisual coronal motions with low signal- 
to-noise ratios of the order of 0. 1%. The technique operates with a set of one-dimensional “surfing” signals extracted 
from position-time plots of several AIA channels through a modified version of Radon transform. The signals are 
used to evaluate a two-dimensional power spectral density distribution in the frequency-velocity space that exhibits 
a resonance in the presence of quasi-periodic PDs. By applying this analysis to the same fan loop structures observed 
in several AIA channels, we found that the traveling velocity of PDs increases with the temperature of the coronal 
plasma following the square-root dependence predicted for slow mode magneto-acoustic waves which seem to be 
the dominating wave mode in the loop structures studied. This result extends recent observations by Kiddie et al. 
to a more general class of fan loop system not associated with sunspots and demonstrating consistent slow mode 
activity in up to four AIA channels. 
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1. INTRODUCTION 

Propagating disturbances (PDs) are systematic translational 
motions of spatially localized emission enhancements in solar 
images. Accurate identification of propagation speeds of PDs 
traveling along coronal loops is of critical importance for clari- 
fying the underlying physical mechanism of this phenomenon. 

Ofman et al. (1997) and DeForest & Gurman (1998), us- 
ing data from the SOHO/ UVCS and SOHO/EYY instruments, 
respectively, reported the first observations of intensity pertur- 
bations along coronal plumes consistent with slow magneto- 
acoustic waves. Signatures of slow mode waves were later de- 
tected in both closed and open loop geometries using a variety of 
instruments, data analysis tools, and simulation techniques (De 
Moortel et al. 2000; Nakariakov et al. 2000; Ofman et al. 1999, 
2000; Marsh et al. 2009; Prasad et al. 2012; Kiddie et al. 2012). 
The onset of the Hinode/ EIS era made the topic of PDs more 
controversial. Some authors have found a clear in-phase relation 
between Doppler shift and intensity oscillations which could be 
attributed to a slow wave propagating upward from the transition 
region into the corona (Wang et al. 2009b, 2009a). Others have 
reported time-variable blueward asymmetry in spectroscopic 
EIS observations, suggesting that the observed PDs might be 
caused by quasi-periodic plasma upflows; therefore, the wave 
theory is not unique (De Pontieu & McIntosh 2010; Tian et al. 
201 la, 201 lb; Warren et al. 201 1 ; Ugarte-Urra & Warren 2011; 
McIntosh et al. 2012; Su et al. 2012). Some of these upflows may 
be coronal counterparts of type II spicules ejecting hot material 
from the chromosphere (De Pontieu et al. 2007, 2011). 

Flows and waves are not necessarily mutually exclusive. 
Quite often, they can coexist inside the same coronal loop, with 
upward plasma flows at footpoints driving slow wave activ- 
ity at higher altitudes (Nishizuka & Hara 2011; Ofman et al. 
2012). In some cases, upward plasma flows are accompanied 


by downward flows reflecting cooling plasma dynamics (Kamio 
et al. 2011; Young et al. 2012; McIntosh et al. 2012), mak- 
ing observational classification of PD signatures even more 
complicated. 

Analysis of the temperature dependence of PD velocities adds 
an important piece of information to this picture by offering 
a potentially crucial test for the slow mode waves in coronal 
loops. If the PDs are caused by a slow magneto-acoustic 
wave, as opposed to a bulk plasma motion, their propagation 
velocities should be controlled by the local sound speed which 
is proportional to the square root of the plasma temperature. 
This dependence can be verified by comparing PD speeds in 
several bands of the extreme ultraviolet coronal emission. Since 
the relative amplitude of PDs tends to be low, determination of 
their speeds over a range of temperatures relies on data analysis 
algorithms enhancing the signal-to-noise ratios (S/Ns). 

Several types of techniques have been developed to reach this 
goal. The initial signal enhancements are commonly performed 
using running-difference image sequences converted into time- 
differenced position-time plots (De Moortel et al. 2000). Ap- 
plying a moving average over a few consequent image frames 
helps achieve an acceptable S /N. Some of the spatial resolution 
is usually sacrificed to obtain a satisfactory confidence level. 
The velocity of the PD fronts can be roughly estimated us- 
ing a simple visual inspection of the position-time plots or by 
applying combined semi-automatic techniques such as identi- 
fying the timings of peaks at each spatial position along the PD 
front and fitting them using a linear regression model (Prasad 
et al. 2012). Any method relying on manual PD processing in- 
troduces a certain degree of subjectivity into the results. More 
elaborate techniques involve a correlation analysis of PD sig- 
natures conducted independently of the observer. By measuring 
the cross-correlation functions between PD signals at each spa- 
tial position, it is possible to evaluate the velocity from the time 
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lag ensuring the highest correlation for a given spatial separation 
(Kiddie et al. 2012; Tian et al. 201 lb; McIntosh et al. 2012). The 
apparent PD speeds can also be estimated by approximating the 
position-time plot with a propagating harmonic wave function 
using a variety of fitting techniques as described by Yuan & 
Nakariakov (2012). A single predefined propagating function 
described by constant values of phase speed, frequency, and 
phase shift limits the range of applicability of this approach. 

To enable a more accurate automated investigation of 
temperature-dependent wave fronts, this paper presents a new 
method for quantifying quasi-periodic PDs in high-resolution 
solar images based on the so-called surfing signals (Uritsky 
et al. 2009) which are obtained by integrating a position-time 
plot along a varying spatiotemporal direction spanning a contin- 
uum of possible propagation speeds. This calculation, which we 
refer to as the surfing transform (ST), leads to an efficient self- 
averaging of random fluctuations preserving the initial time res- 
olution of the image sequence without distorting the frequency 
structure of the PDs. A two-dimensional (2D) Fourier power 
spectrum of the ST makes it possible to detect low-amplitude 
PD events as double resonances in the frequency-velocity space. 
We derive an analytical fit for the velocity resonance which can 
be used for measuring the scale length, the amplitude, and the 
propagation speed of the disturbances. 

The developed technique makes it possible to quantify 
counter- streaming signals traveling simultaneously in opposite 
directions within the same loop structure. Our method does not 
rely on a particular wave form of the PD signal and can be ap- 
plied to both harmonic and nonlinear propagating features. It can 
also be applied to repetitive propagating oscillations covering a 
wide range of frequencies without a well-defined characteristic 
periodicity. The primary limitation of the method is that it loses 
spatial information of the disturbance converting a spatiotem- 
poral PD pattern into one-dimensional temporal signals. In its 
present version, the transform is intended to investigate average 
PD characteristics and their time variability, but not detailed spa- 
tial inhomogeneities associated with propagation besides first- 
order damping effects incorporated into the derived fit. 

The ST algorithm is applied to multi-temperature coronal im- 
ages obtained from the Atmospheric Imaging Assembly (AIA) 
on board the Solar Dynamics Observatory (SDO). We study the 
fan loop structure of the non-flaring solar active region NOAA 
AR 11082 which exhibits propagating emission oscillation in 
several AIA channels, some of them with rather high noise lev- 
els. The results show a clear square-root dependence of the PD 
propagation speed on the loop temperature, consistent with the 
slow wave model. 


2. METHODS AND DATA 

2.7. Surfing Transform Technique 

We define the ST as a time-dependent mean a sit, v) of 
the position-time plot a(x , t) computed along the straight line 
running through the point x = 0 with the slope v as shown in 
Figure 1: 


V rt+L/v 

as(t,v)=— / a(x = xs, t')dt ' . (1) 

^ Jt 

Here, t is the starting time, xs = v x (F — t) is the averaging 
path, L is the propagated distance, and v is the surfing velocity. 
The calculation of the ST is conceptually close to the 2D Radon 
and Hough transforms (see, e.g., Jones & Davila 2009, and 
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Figure 1. Schematic drawing illustrating the extraction of the ST signal from a 
position-time plot. Dashed lines show the fronts of a periodic PD described by 
phase speed vo. The solid line is the average path defined by the surfing speed 
v and the starting time t. L is the size of the region studied. 
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references therein); however, its application is different. Instead 
of using the integral transform (1) for enhancing position-time 
plots, we focus on its resonance behavior. 

If a(x , t) is a periodic disturbance propagating with speed vo, 
the dynamic range of as maximizes at v = vq. Consider the ST 
of a one-dimensional harmonic wave oscillation: 

v rt+L/v 

as(t, v) = — ao cos(kov x (F — t) — coot')dt' 

L Jt 

a^v 

= — -{sin([k 0 f - co 0 ][t + L/v] - k 0 vt) 

L(k 0 v - (o 0) 

— sin([koP — coo]t — kovt)} 

sin(§) / , , .x ... L(k 0 v-co 0 ) 

= ao cos(cD()t + §) with § = , 

§ 2v 

( 2 ) 


in which we used the identity sin(a) — sin(/3) = 

2 8 ^( 0 ? — /?/ 2 )cos(a + jft/ 2 ). In the dispersionless case, the 
phase speed vo = coo/ko and, therefore, the velocity detun- 
ing parameter § = Lcooiv — vq)/(2vvo). As v approaches vo, § 
tends toward zero causing a sit, v) to converge to aix = 0, t). 
In this limit, the averaging paths in (1) become constant phase 
lines yielding the largest possible ST amplitude given by ao. 

Since (1) is a linear procedure, the above result also applies 
to a sum of multiple harmonic modes. To investigate more 
complex processes, we calculate the ST spectrum describing 
the 2D distribution of the power spectral density over a range of 
Fourier frequencies and surfing velocities: 

v I rt+L/v _ 2 

Psiv,f) = — \ I / aiv x it'— t), t')e~ l2n ^ dt'dt , 

L 17-00 Jt 

( 3 ) 

as well as the velocity spectrogram P s iv, f = /o) for a 
fixed frequency /q. The shape of the velocity spectrogram of 
a harmonic wave at /o = coo/i2n) is given by the square of 
the sine (§) function (see Equation (2)) reaching its maximum at 
v = v 0 : 

P s = Po sinc 2 (L&>oO - vo)/(2vv 0 )), (4) 

where Po is the height of the spectral peak. As the propagation 
distance L and the frequency coo of the oscillation increase, 
the peak becomes narrower and the resonance measurements 
become more precise. The increase of vo has the opposite effect, 
leading to a broader Psiv) peak. 

2.2. Statistical Testing 

Panels (a)-(e) of Figure 2 show the application of the 
described technique to three superposed wave signals mixed 
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(a) 




Figure 2. Testing the ST technique on upward and downward PDs. (a) 
Position-time plot of the simulated wave signal defined by Equation (5). 
(b) Stack plot of the signal at five selected x positions, (c) Surfing signals 
obtained using three values of v matching the phase speeds of the wave 
components, (d) ST spectrum computed using 400 ST signals covering the 
velocity range v e [—200, 200] km s -1 . (e) Velocity spectrograms representing 
horizontal sections of the ST spectrum at /o = 5.00 mHz (dotted blue line) and 
/o = 3.3 mHz (solid black line). Note that the method recognizes concurrent 
downward and upward motions in the same or different frequency bands. 

(A color version of this figure is available in the online journal.) 


signal with negative velocity propagates in the negative v di- 
rection. The surfing signals in panel (c) are extracted from the 
position-time plot (a) by tuning v to the three wave speeds. 
The ST spectrum shown in panel (d) features three peaks cor- 
responding to these waves. The velocity spectrograms repre- 
senting horizontal cuts of the ST spectrum through the wave 
frequencies (two of which coincide) are provided on panel (e). 
For a more precise result, the spectrograms were averaged over 
the frequency bands shown in panel (d) by dashed horizontal 
lines. The spectrograms demonstrate surfing resonances at the 
expected v values. The negative velocity peak is easily identifi- 
able, demonstrating the ability of the method to recognize PDs 
traveling downward and upward along coronal loops during the 
same observation interval. This capability constitutes one of the 
main advantages of the developed method compared to tech- 
niques relying on cross-correlation analysis (Kiddie et al. 2012) 
or nonlinear regression fits (Yuan & Nakariakov 2012) using 
single-component wave models. 

The phase speed of the slow mode waves in coronal loops 
is controlled by the local sound speed c 2 s = yp/ p oc T, where 
y is the adiabatic index, p and p are the kinetic pressure and 
the mass density of the coronal plasma, and T is the plasma 
temperature. Assuming that the wave propagates in a bundle of 
isothermal coronal loops excited by the same periodic process 
in the transition region but described by different values of T, 
we expect the PD velocities in these loops to scale as the square 
root of T while the frequencies remain approximately the same. 

The velocity-temperature scaling vq a T 1/2 , indicative 
of the slow mode, can be tested based on ST analysis of 
different SDO AIA channels with distinct temperature peaks 
of sensitivity (Boerner et al. 2012; Lemen et al. 2012). Figure 3 
presents a numerical simulation of such a multi-temperature 
measurement in the presence of additive Brownian noise. The 
first column shows position-time plots of harmonic wave 
oscillations modeling slow waves with vq = 90 km s -1 
observed in the 171 A channel. The noise intensity increases 
from top to bottom. Similar a(x , t) plots are constructed for 
the 131 A, 193 A, and 211 A AIA channels, with vo satisfying 
the temperature scaling of the slow magneto-acoustic wave 
(68 km s -1 , 136 km s -1 , and 144 km s -1 , respectively). The 
second column shows the velocity spectrograms P$( v) of the 
channels computed at the base wave frequency / 0 = 3.0 mHz. 
Parameters vo, L , and Pq were adjusted to obtain the best 
least-square fits (4) plotted with solid black lines. Solid and 
dashed vertical lines show the estimated and true vo values, 
respectively. It can be seen that our method yields reasonable 
velocity estimates for quite low S/Ns. 

The last two rows of panels in Figure 3 show the results 
of a null testing of our method using Brownian and white noise 
models without admixtures of periodic signals. In both cases, the 
ST spectrograms reveal no well-defined velocity peaks observed 
in the presence of PDs. 

2.3. SDO AIA Observations 


with a correlated Brownian noise (Mandelbrot & VanNess 1968) 
resulting in a 10% S/N: 

3 

a(x, t) = E oq cos {kiX -(0it) + a$ Xtt , (5) 

i=i 

in which [oq, <z> 2 , &>3]/(27 t) = [5.00,3.33,3.33] mHz, 

[v\, V2, U 3 ] = [—100, 50, 120] km s -1 , and ki = cc>i/Vi. The 


The described methodology was applied to SDO AIA images 
of solar active region NOAA AR 11082 collected on 2010 
June 19 during 4:00-10:00 UT with an angular resolution of 
~ 0.6 arcsec and a maximum available cadence time of ~12 s. 
The region was located in the northern hemisphere close to 
the disk center. It had a bipolar structure not associated with 
a sunspot, and showed no flaring activity above B -class for 
the time period studied. We use level 1.5 data from the cutout 
service. The images were derotated by the SHIFT_IMG routine 
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Figure 3. ST analysis of simulated wave signals in four SDO AIA channels, as compared to non-propagating random noises. Left: time-differenced position-time 
plots of the 171 A channel. Right: velocity spectrograms obtained for each channel; the phase speed of the wave is adjusted to the temperature of the channel assuming 
that the wave is slow magnetosonic. Color coding: 131 A — blue, 171 A — red, 193 A — green, 211 A — brown. Black curves are least-squares fits using the theoretical 
model (4), solid (dashed) vertical lines are the measured (predicted) propagation velocities for each channel. The amount of additive noise increases from top to bottom 
as shown by the provided power signal-to-noise ratios (S/N). The tilted red lines on the left panels show the phase speeds and wave periods measured at 171 A. The last 
two rows of panels corresponding to the cases of random fluctuations with Brownian and “white” spatial spectra showing no wave PD signatures on ST spectrograms, 
as expected. 

(A color version of this figure is available in the online journal.) 


of the SolarSoft package, and co-aligned as described by Viall & 
Klimchuk (2011, 2012). As can be seen from Figure 4, NOAA 
AR 1 1082 has an extended system of fan loops. The PD activity 
in these loops was studied using the 131 A, 171 A, 193 A, 
and 211 A channels. The remaining AIA channels showed 
insufficient count rates and were excluded from the analysis. 

Three of the studied AIA passbands (131 A, 193 A, and 
211 A) are known to be described by multithermal response 
functions (Del Zanna 2012; Del Zanna & Storey 2012; Del 
Zanna et al. 2012) which require special attention. The 131 A 
channel reflecting the contribution from Fevm, Fexx, and 
Fe xxiii lines has two distinct peaks of sensitivity at about 
0.56 MK and 10.8 MK. The high-temperature peak is ~2.5 times 
lower than the first peak (Boerner et al. 2012; Lemen et al. 
2012). The significance of the high-temperature 131 A peak for 
our analysis is limited for two reasons. First, the studied non- 
flaring coronal region is relatively cool and most likely contains 
no ~10 7 K plasma. Second, even if such plasma exists, the 
expected slow mode speed at this temperature is beyond the 


range of measurable speeds. According to DeForest & Gurman 
(1998) and subsequent studies (see, e.g., Aschwanden 2006 for a 
review), the characteristic traveling velocity of acoustic waves in 
coronal structures observed in the approximately monothermal 
171 A channel ( T ~ 0.9 MK) is of the order 100 km s -1 which 
translates into ~350 km s -1 at T = 10.8 MK. On the 24 s 
Nyquist time scale (two AIA sampling intervals), such a wave 
would travel a distance of ~8.5 Mm, which is comparable to 
the projected length of the studied loop region and cannot be 
accurately detected by any method. Based on these arguments, 
we consider the high-temperature peak of the 131 A emission 
to be irrelevant to our analysis and focus on the main sensitivity 
peak of this channel. 

The 193 A passband is also multithermal, reflecting a strong 
contribution from Fevm, Feix, Fexi, and Fevn lines, as 
well as a much weaker emission from lower-temperature lines 
dominated by O V and Fe vn (Del Zanna et al. 2011; Kiddie et al. 
2012). The main cool contribution to the 193 A passband comes 
from Fevm and Feix lines. As shown by Kiddie et al. (2012), 
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Figure 4. (a) Snapshot of the active region NOAA AR 1 1082 as seen in the 171 A SDO AIA channel at 7:00:00 UT. The dashed white polygon shows the domain of the 
PD activity studied in Figures 5 and 6; the solid white line outlines the adjusted domain used to analyze two selected PDs reported in Table 1. (b) Multi- spectral AIA 
images of the same region obtained around the same time; temperature values are for the peaks of the response functions based on the 2013 release of the CHIANTI 
atomic database (Dere et al. 1997; Landi et al. 2013). 

(A color version of this figure is available in the online journal.) 


it can reach up to 40% of the total emission in this channel in 
sunspot regions and in that case cannot be ignored. According 
to the same study, the significance of the cool 193 A component 
is much lower in non-sunspot regions such as NOAA AR 1 1082 
where it can be responsible for about 15%-25% of the observed 
emission. Due to the smallness of this contribution, we do not 
address it directly, keeping in mind that the estimated speeds 
can be somewhat slower than what might be expected for the 
main 193 A sensitivity peak centered at 1.6 MK. As shown in 
the next section, such bias does exist and it is small compared 
to the observed temperature dependence. 

The low temperature emission of the 21 1 A passband is gener- 
ally a factor of two smaller than the corresponding contribution 
of the 193 A passband (Boerner et al. 2012). Taking o this fact 
into account, we associate the temperature of the 21 1 A channel 
with its main response function peak around 1.8 MK. 

3. ST ANALYSIS RESULTS 
3.1. PD Activity in the Selected Fan Loop Region 

We found that many of the studied fan loops sectors carried 
quasi-periodic PDs during the entire 6 hr interval of obser- 
vations. The disturbances travel predominantly in the upward 
direction away from the base of the corona, with a periodicity 
of 4-8 minutes and an apparent propagation speed of ~40 to 
~180 km s _1 . Here we focus on the fan loop sector highlighted 
by dashed white lines in Figure 4(a) which shows the most stable 
PD activity. 

Figure 5 shows the results of the ST analysis of the chosen 
sector. The position-time plots shown in the left panels exhibit 
PD patterns. To enhance the PD signals, the plots were subjected 
to temporal differencing with a 60 s time lag. The pairs of 
straight lines added to each panel show measured PD velocities, 
with the interval between the lines showing the characteristic 
relaxation time / 0 _1 ~ 310 s of PD fronts given by the highest 
resolvable harmonic of the ST spectrum in the least noisy 17 1 A 


channel. The ST spectra of 171 A emission typically contain 
several harmonics caused by the intermittent structure of the 
studied signals. The highest harmonic used here represents the 
characteristic duration of single intensity perturbations and is 
the most relevant to the velocity analysis of the individual PD 
fronts. 

The velocity spectrograms fitted with Equation (4) are pre- 
sented in the central panels. The obtained velocity estimates 
(marked by solid vertical lines) are approximately consistent 
with the slow mode speeds (dashed lines) predicted based on 
the measured 171 A channel speed. The resonance peaks are 
broader for larger vq values, in agreement with the analytical 
dependence (4). 

The estimated scale length L varies between ~ 1 1 Mm at 1 3 1 A 
and 20-21 Mm at 171 A and 193 A wavelengths. This behavior 
could be a mixture of several factors, such as e.g., different levels 
of noise in the AIA channels, the temperature-dependent scale 
height of sound waves in a radially diverging coronal magnetic 
field (e.g., Torkelsson & Boynton 1998), and viscous dissipation 
of wave energy leading to pronounced altitudinal stratification 
of the slow mode amplitude (Ofman et al. 2000). 

The right column of panels in Figure 5 shows surfing signals 
obtained from the respective a(x , t) plots using the estimated 
PD speeds. The shape of the signals confirms the periodic na- 
ture of the observed activity. The apparent cross-channel coher- 
ence of the signals, which represents temporal oscillations at 
the base (x = 0) of the studied wave region, suggests that the 
wave fronts in the studied multi-temperature loops were excited 
nearly simultaneously at the base of the corona, possibly by 
a global p - mode oscillation. The coherence is the most pro- 
nounced for the strongest wave fronts. For instance, the two 
intense wave fronts marked by dashed vertical lines are seen 
at approximately the same times in all AIA channels. We con- 
ducted a separate analysis of this pair of wave fronts within a 
narrow loop sector (solid white boundary in Figure 4(a)) where 
the fronts were the most intense. The position-time plots of the 
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Figure 5. Analysis of quasi-periodic PDs observed in the fan loop sector outlined by the dashed white line in Figure 4(a): position-time plots (left), ST-based velocity 
spectrograms (center), and ST signals constructed using vo values measured for each passband. Black solid curves overplotted with velocity spectrograms are the best 
least- squares fits based on Equation (4) yielding the estimates of the propagation velocity vo and the projected PD scale length L. Dotted lines show the shape of each 
spectrogram under the assumption of temperature-independent wave propagation at the 171 A phase speed. Solid vertical lines show the measured PD speeds; dashed 
vertical lines are the predicted slow magnetosonic wave velocities based on the 171 A measurement. Vertical lines on the ST signal mark the timing of the two selected 
fronts in the 171 A channel discussed in the text. The noisiest 131 A channel is represented by the original (gray) and the low-pass filtered (black) ST signals. 

(A color version of this figure is available in the online journal.) 


Table 1 

Velocity Ratios Measured for the Entire Wave Event and Two Selected Fronts, 
as Compared to the Predicted Sound Speed Ratios 


Wavelength 

All Fronts 

Selected Fronts 

Predicted 

131 A 

0.63 ± 0.03 

0.86 ± 0.04 

0.79 

171 A 

1 

1 

1 

193 A 

1.23 ±0.08 

1.31 ±0.06 

1.34 

211 A 

1.50 ±0.23 

1.50 ±0.15 

1.42 


Note. See text for details. 


fronts were replicated several times to enable the application of 
our ST algorithms to episodic PD activity. 

Table 1 reports the velocities of the selected replicated fronts 
normalized by the velocity in the 171 A channel, along with 
the normalized velocities obtained for the entire time interval 
shown in Figure 5. The errors are the propagated standard 
deviations from the nonlinear least-squares fit (4) with the 
Gaussian weighting. For reference, the last column of the table 
provides the velocity ratios for the sound wave. Comparing 
velocity ratios has the advantage of eliminating the uncertainty 
associated with projection effects and the unknown value of y 
which can vary between 1 and 5/3 (Van Doorsselaere et al. 
2011). It can be seen that the velocity ratios measured by our 
method are fairly close to the prediction for the slow mode 
wave. The agreement with the theory is somewhat better for the 
selected intense fronts due to their higher S /N. 



-200 -100 0 100 200 
v, km/s 


Figure 6. Bidirectional velocity spectrograms obtained by integrating ST spectra 
over the frequency range 0.002-0.02 Hz (periods 50-500 s, or 0.83-8.3 minutes). 
The 171, 193, and 211 A passbands show a prevalence of upwardly moving 
features. The 131 A spectrogram is approximately symmetric due to low count 
rates in this AIA channel as can be seen from the inset plot displaying non- 
normalized Ps(v) plots. 

(A color version of this figure is available in the online journal.) 

Our analysis revealed no signatures of downwardly propa- 
gating PDs in the studied loop coexisting with temperature- 
dependent upward motion. Figure 6 presents bidirectional veloc- 
ity spectrograms which were integrated over a frequency range 
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Figure 7. Results of the PD survey using the ST technique, (a) 171 A image divided into sectors. The color-coded bars show the presence of identifiable PD activity 
in the corresponding AIA channel, (b) Velocity ratios between the measured PD velocities and the velocity in the 171 A passband vs. loop temperature given by the 
maximum SDO AIA sensitivity in each channel. The theoretical square-root dependence expected for slow mode waves is plotted by the black solid line; symbols 
show the results for the fan loop sectors indexed in panel (a) in which PDs were observed at one or more AIA wavelengths including 171 A. The dashed line represents 
the average velocity ratio characterizing different fan sectors, with vertical error bars showing plus minus one standard error. 

(A color version of this figure is available in the online journal.) 


(b) 



(a) 



0.002-0.020 Hz to take into consideration broadband emission 
oscillations. This range includes the PD activity described above 
as well as other possible motions associated with flows and 
waves. It can be seen that the spectral power in the 171, 193, 
and 211 A passbands associated with the negative propagation 
direction is a factor of 4-5 smaller compared to the power in the 
positive (upward) direction and shows no distinct maximum at 
Vo < 0. The much lower count rates of the 131 A channel leave 
a possibility of an unresolved downward motion in this channel. 
Overall, Figure 6 suggests that even if an undetected downward 
motion is present its intensity is negligibly small compared to 
that of upwardly traveling disturbances. 

3.2. Distribution of PDs across Other Fan Loops 

Using the same methodology, we studied the spatial distri- 
bution of PD activity across the loops belonging to the same 
active region. The fan loop system of NOAA AR 11082 was 
divided into 50 narrow sectors diverging radially from the cen- 
ter of the fan loop structure and characterized by equal angular 
(~2?3) and radial (~17.3 Mm) sizes. The sectors were approx- 
imately parallel to the local magnetic field traced by the fan 
loops. For each sector and the AIA channel, we constructed a 
position-time plot describing the spatial distribution of the AIA 
intensity along the sector as a function of time during the same 
interval as that investigated in Figure 5. To enhance periodic 
signals, the position-time plots were subject to 60 s temporal 
differencing as before. 

The results of the survey suggest that upward PD activity 
is ubiquitously present throughout the studied fan system as 
shown in Figure 7(a). The results of the velocity measurements 
in different sectors are provided in Table 2. These measurements 
confirm the tendency of the PD velocity to increase with the 
temperature of the coronal plasma. Most of the sectors exhibited 
sufficiently strong PDs in at least two AIA channels making it 
possible to compute the velocity ratios in the individual fan loop 
sectors and compare the estimated values with the slow wave 
mode prediction. Panel (b) of Figure 7 presents the results of 
this verification. Here, the velocity ratios estimated for several 
loop sectors are plotted using various symbols as shown in 
the legend. The plotted values were aggregated into groups 


Table 2 

Parameters of PD Activity in the Fan Loop Sectors Shown in Figure 7 


Slice Index 

f, mHz 


vo, km s 1 




131 A 

171 A 

193 A 

211 A 

0-4 






5-9 

1.42-3.19 

47 ±3 

71 ±5 

83 ± 13 

95 ±4 

10-14 

2.42-5.43 


95 ±4 

112 ±6 


15-19 

2.70-4.56 

49 ±6 

82 ±3 

132 ± 17 

185 ±98 

20-24 

1.94-4.36 



115 ± 14 


25-29 

1.66-3.73 


100 ± 10 

122 ± 14 

98 ±7 

30-34 

1.83-2.63 



116 ± 10 

165 ± 57 

35-39 

40-44 

2.01-2.89 

64 ±3 

86 ±3 



45-49 

3.02-4.35 


72 ±4 

95 ±7 



of five adjacent sectors to improve statistical accuracy. The 
dash-dotted line represents the average between these values 
for each temperature. 

The observed temperature scaling is in a reasonable agree- 
ment with the theoretical square-root dependence plotted with 
the solid black line. This tendency is statistically significant for 
the 193 A and 131 A velocity ratios but not for the 21 1 A ratio 
due to larger propagation velocities observed in this channel 
(Table 2) yielding a broader Ps peak (see Equation (4)) leading 
to a greater statistical uncertainty. It is interesting to note that 
the 193 A velocity ratios are systematically below the theoret- 
ical slow mode dependence. This bias likely reflects the cool 
emission component present in this passband. The 131 A ve- 
locity ratios are also systematically smaller than the predicted 
values, but this is most probably due to the high level of noise in 
this channel leading to false Ps peaks at vo ~ 40-50 km s -1 as 
suggested by statistical tests presented in Figure 4. Such peaks 
reflecting the correlation structure of non-propagating fluctua- 
tions can interfere with the actual velocity resonances causing 
an asymmetric broadening of Ps peaks toward lower PD speeds 
and yielding underestimated Vo values. Some of the bias can 
also be caused by the contribution from the cool branch of 
the 131 A temperature response (Tian et al. 2011b; McIntosh 
et al. 2012). The magnitude of these effects, which could in 
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principle be removed using a more sophisticated fitting func- 
tion, is nevertheless smaller than the observed square-root 
temperature trend. 

3.3. Relationship with Flows 

Our results strongly suggest that the studied fan loop system 
is dominated by upwardly moving disturbances which are 
consistent with the slow magneto-acoustic wave model. The 
unavailability of high-cadence spectroscopic observations for 
these events makes it impossible to verify this interpretation 
using Doppler measurements as was done by other authors, 
see, e.g., Wang et al. (2009b, 2009a). However, the evidence 
obtained presents an explicit case in favor of slow mode waves 
in the studied solar region. 

It should be clarified that the observed temperature depen- 
dence of PD velocity cannot be associated with a superposition 
of downward and upward bulk motions considered in some pre- 
vious works. A number of recent studies reported unambiguous 
signatures of such flows using conjugate coronal imaging and 
spectroscopic observations. Some of these works were already 
mentioned in Section 1. In particular, De Pontieu & McIntosh 
(2010) identified ubiquitous and faint upflows with characteris- 
tic speeds 50-150 km s -1 in coronal loops associated with plage 
regions. These flows caused blueward asymmetries of spectral 
line profiles in footpoint regions of coronal loops and could 
be related to high-temperature plasma jets originating in the 
chromosphere. The hot plasma upflows can be accompanied by 
downflows of radiatively cooling material closing the mass cycle 
of the coupled chromosphere — corona system (McIntosh et al. 
2012). The presence of counterstreaming flows complicates the 
interpretation of line profiles that can exhibit a net Doppler shift 
asymmetry suggestive of an apparent temperature-dependent 
plasma speed. 

We note, however, that such contraflows would be resolved 
by the ST technique if they were present in the fan loop system 
studied here. As shown in Section 2, the ST analysis is able 
to differentiate between counter-PDs and would selectively 
detect both types of PDs if they were simultaneously present at 
either different or the same coronal altitudes. However, no such 
periodic contraflow events (McIntosh et al. 2012) were found in 
the fan loop regions studied. The 1 3 1 A passband, which exhibits 
a roughly symmetric ST velocity spectrogram (Figure 8), has 
very low count rates and is likely to be dominated by noise. In 
some locations, we did observe occasional features descending 
toward the transition region such as the one shown in Figure 8. 
A closer inspection has shown that these features are caused by 
proper motion of several closed loop filaments relative to the 
steady fan loop system. This effect is similar to the intensity 
oscillations associated with Alfvenic waves observed by Tian 
et al. (2012), which cause loop displacement rather than density 
variations, and it cannot be attributed to plasma downflows. 

Obviously, the existence of plasma flows cannot be ruled out 
based on the imaging diagnostics alone. Based on our measure- 
ments, it is unlikely that bulk plasma flows play a leading role 
at the coronal altitudes covered by our analysis. The observed 
emission dynamics indicate that (1) the motion occurs in only 
one direction and (2) it exhibits wave-like temperature scaling. 
However, there is a possibility that significant flow motions un- 
resolved by our method exist near the footpoints of the fan loops 
studied. These motions could be synchronized by p - mode os- 
cillations and be part of the global chromosphere-corona mass 
cycle as proposed by McIntosh et al. (2012). At the same time, 
they could excite slow mode and other magnetohydrodynamic 



time, UT 

Figure 8. Example of apparent downward motion in the lower portion of the 
fan loop system studied coexisting with upward PDs at higher altitudes. Panel 
(a) shows the analyzed region outlined by a white dashed polygon; panel (b) is 
the detrended position-time plot representing this region. The downward PDs 
in this event are an artifact caused by the displacement of a set of descending 
closed loops marked with a dashed line on panel (a) which cross the fan loop 
filaments included in the virtual slit. 

(A color version of this figure is available in the online journal.) 


waves in the upper coronal regions. In this scenario, the faster 
propagating upward flows would provide energy for the waves. 
In open field line geometry, waves can carry a sizable portion of 
this energy away from the transition region into the solar wind. 
A similar flow-wave interaction has been recently reproduced 
in magnetohydrodynamic simulations of a coronal active region 
(Ofman et al. 2012; Wang et al. 2013). More experimental and 
theoretical studies will be needed to verify the validity of the 
scenario and its underlying physics. 

4. CONCLUSIONS 

We have presented a new methodology for measuring pa- 
rameters of PDs in multi-temperature coronal images. The key 
element of the proposed approach is the analysis of surfing sig- 
nals which demonstrate a resonance behavior in the presence 
of quasi-periodic PDs. The developed methodology has been 
applied to the analysis of PDs in coronal fan loops in AR 1 1082 
across a range of temperatures covered by four SDO AIA 
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channels. We found that the traveling velocity of PDs propagat- 
ing along the loops obeys a square-root temperature dependence 
predicted for slow mode magneto-acoustic waves which seem 
to be the dominating wave mode in the loop structures studied. 

Recently, Kiddie et al. (2012) arrived at a similar conclusion 
based on a detailed analysis of PDs in coronal fan loops 
anchored at sunspots. Using three AIA channels, they found 
that the velocities of such PDs depend on the temperature as 
predicted for the slow wave mode. However, they reported no 
conclusive observations for PDs in non-sunspot locations. Our 
study extends their findings by showing that the temperature- 
dependent slow mode activity can be observed simultaneously 
in four AIA channels in a non-sunspot active region. It remains 
to be clarified whether or not this wave activity is driven by 
undetected fast plasma flows at the base of the corona as has 
been proposed by some authors (De Pontieu & McIntosh 2010; 
Tian et al. 2011a, 2011b; Warren et al. 2011; Ugarte-Urra & 
Warren 2011; McIntosh et al. 2012; Su et al. 2012; Ofman et al. 
2012 ). 

The range of applicability of the technique presented here 
involves a variety of more general tasks beyond the scope of 
this paper. For instance, it would be of interest to apply the 
technique to plume regions. Plumes were the first coronal struc- 
tures where slow magnetosonic waves were identified (e.g., Of- 
man et al. 1997, 1999; DeForest & Gurman 1998). Recently, 
Tian et al. (2011c) found no clear temperature dependence of 
PDs observed in plume-like structures rooted in various solar 
regions including polar and equatorial coronal holes as well 
as the quiet Sun. The velocity measurements conducted in 
their study were based on visual inspection of position-time 
plots. The technique described in our paper can be used to 
obtain more accurate PD velocity estimates in such structures 
and to test earlier conclusions regarding the presence or ab- 
sence of the magneto-acoustic v — T scaling in multithermal 
plumes. 

In principle, ST analysis can be applied to any superpositions 
of multiple PD processes characterized by different velocities 
and/or frequencies as illustrated in Figure 2, including strongly 
nonlinear wave forms. It could be especially helpful for ana- 
lyzing wave damping effects since the ST algorithm provides 
self-consistent evaluation of the amplitude and the scale length 
of the wave. Such measurements can be conducted in various 
astrophysical contexts not limited to the solar corona. 

The authors thank James Klimchuk and Tongjian Wang 
for useful discussions. V.U. was supported by the NASA 
Grant NNG11PL10A 670.002 through the CUAs Institute for 
Astrophysics and Computational Sciences. L.O. was supported 
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